MZ-TH/95-21 
TECHNION-PH 97-02 



Inclusive Semileptonic Decays of Polarized A5 
Baryons into Polarized r-Leptons 



ON 

: S. Balk*, J.G. Korner 1 !] and D. Pirjol s 

a 1 



X 



1 Johannes Gutenberg- Universitat, Institut fiir Physik, Staudingerweg 7, 

D-55099 Mainz, Germany 



■ Department of Physics, Technion-Israel Institute of Technology, Haifa 32000, Israel 

m 
m 
o 
t> 

ON 

a,: 



Abstract 

We employ OPE techniques within HQET to calculate the inclusive 
semileptonic decays of polarized A& baryons. Lepton mass effects are in- 
cluded which enables us to also discuss rates into polarized r-leptons. We 
present explicit results for the longitudinal polarization of the r in the A& 
rest frame as well as in the (r~,u T ) cm. frame. In both the A& rest 
frame and in the (t~ ,P T ) cm. frame we make use of novel calculational 
techniques which considerably simplify the calculations. The transverse po- 
larization components of the r are calculated in the (t~,u t ) cm. frame. 
We delineate how to measure the full set of 14 polarized and unpolarized 
structure functions of the decay process by angular correlation measure- 
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ments. A set of observables are identified that allow one to isolate the con- 
tributions of the two 0(l/m|) nonperturbative matrix elements Kb and e&. 
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I. INTRODUCTION 



Large samples of the bottom baryon have been produced on the Zq at LEP and 
are expected to be produced in future colliders. Advances in microvertexing techniques 
have allowed for efficient means of A& identification. For example, when LEP was running 
on the Z one had ~ 2.2 x 10 5 bb pairs per 10 6 Z-decays. Of these approximately 10% go 
into A;, baryons of which again « 20% decay semileptonically (e, \i and r). Thus one can 
expect a sample of 4000 inclusive semileptonic A& (or Aj,) decays for every 10 6 Z-decays. 
Plans at the SLC call for altogether 3 x 10 6 produced Z's. The quality of the A& data 
from the SLC will even be better because its small beam size provides for an excellent 
definition of the A& production vertex. 

The A fe 's produced on the Z-peak are expected to be quite strongly polarized Q. 
This calls for a consideration of polarization effects in the A& decays which could be 
used to determine the polarization of the A&. Approximately 10% of the total (e + 
fi + t) semileptonic decay sample have a r-lepton in the final state. The r-lepton is 
sufficiently heavy which necessitates the inclusion of lepton mass effects in the dynamical 
rate calculations, apart from pure phase space effects. We will study also the polarization 
of the r-lepton which, because of mass effects, differs from its naive mi = limiting 
value. The r-polarization can be experimentally determined from its subsequent decay 
distributions. 

In this paper we present all the necessary tools to calculate inclusive semilep- 
tonic decays of polarized A^'s including r-lepton polarization effects and spin-spin, spin- 
momentum and momentum-momentum correlation effects. Our calculation makes use 
of the heavy quark effective theory (HQET) and the operator product expansion (OPE) 
method as applied to heavy hadron decays. There is some overlap of our work with earlier 
results on inclusive A& decays in |^||[|,|5]||[7||| and the more recent paper by M. Gremm, 
G. Kopp and L.M. Sehgal on polarization effects in inclusive semileptonic A^-decays @. 
At the technical level we differ from the analysis of the above papers in that we employ 
helicity techniques to derive compact forms for the differential decay distributions includ- 
ing polarization effects. Also, by using suitable phase space variables, the integration of 
the higher order terms in the operator product expansion become much simplified, i.e. 
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there are no surface term contributions. 

In Sec. II we list our results on the hadronic matrix elements needed for the subse- 
quent semileptonic rate calculations. As the techniques of deriving the hadronic matrix 
elements are quite standard by now we do not dwell much on theoretical background but 
immediately proceed to the final results which we list in terms of a set of spin-dependent 
invariant structure functions. In Sec. Ill the invariant structure functions are related to he- 
licity structure functions which determine the complete angular structure of the polarized 
decay distributions involving the three helicity angles of the decay process. We write down 
the full differential decay distribution of polarized A b -decays into longitudinally polarized 
r-leptons. We give analytical and numerical results on decay distributions integrated over 
phase space. In Sec. IV we discuss the case of transversally polarized r-leptons. 

Sec.V is dedicated to a calculation of the longitudinal polarization of the r in the 
A;, rest frame as well as of its azimuthally averaged transverse polarization component 
in the plane spanned by the r and the polarization vector of the A&. The calculation is 
considerably simplified by extracting the rate as the absorptive part of the appropriate 
one-loop contribution to the fA b — > fA b scattering amplitude. Sec. VI contains our con- 
clusions. More detailed results have been collected in two Appendices. In Appendix A we 
list results on go-integrated helicity structure functions. In Appendix B we give fully in- 
tegrated analytic results for the seven structure functions that describe the angular decay 
structure of unpolarized A b decays. 

II. HADRONIC MATRIX ELEMENTS 

The dynamics of the hadron-side transitions is embodied in the hadronic tensor W^ u 
which is defined as 

W^(q , q 2 ,s) = (2nfY,S\Pi-V-Pz)(MPi,s)\jf\Xj(P^ (1) 

x 

where Jj*(j = c,u) is the hadronic current inducing b — > c and b — > u transitions. The 
hadron tensor is a function of two kinematic variables which we choose as go and g 2 . Note 
that we are not summing over the A^-spin such that the hadron tensor depends also on 
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its spin four-vector s. The structure of the hadron tensor can then be represented by an 
expansion along a standard set of covariants [10| (v = pi/rri\ b ) 

= -g^Wi + v»v v W 2 - ie^v^Wa + q^q u W 4 + [<fv v + v"q v )W 5 



+ 



— q.s 



-gi»G x + v»v v G 2 - ie^v a q p G 3 + <fq u Gt + (q»v u + q v v»)G t 



+ sV)G 6 + {s»q v + sV)G 7 



+ier#v a s f} Gs + ie! tt "#q a s f} G. 
+{s»v v - s v v»)G l0 + (s»q v - sV)Gu 
+( v *(r>&r qaVj3 s + v v e IM ^q a VpS 1 )G 12 



+{q^ a ^q a v pSl + q u t^q a Vp Sl )G^ . 



(2) 



The last four invariants Gio, Gn, G\ 2 and G\-& are so-called T-odd invariants which are 
fed by CP-odd and/or imaginary part contributions. They are zero for Standard Model 
couplings when loop effects are neglected. They will therefore be disregarded in the 
following. Note that the invariants W4, W5, G4, G$, G7, Gn and G13 do not contribute to 
inclusive semileptonic decays in the zero lepton mass case. Since we are also considering 
decays involving the massive r-lepton we must keep the full set of invariants implied by 
Eq.(§. 

In order to compute the spin- independent structure functions Wi,...,W^ and the 
spin-dependent structure functions G\, G 9 we resort to the well-known OPE techniques 
in HQET. The requisite steps in this calculation are so well documented in the literature 
P fI^jn| , p^JT3| , p!4] , p!5| , p!6| , p!7| that we can forgo a description of the intermediate steps and 
immediately list the final result of the OPE analysis. The hadron tensor is obtained as 
the absorptive part of the forward matrix element (W^ = — ^Im T^ v ) 



T^(q ,q 2 ) = -i(A b ( Pl ,s)\ / d 4 xe-^ x TJ^(x).F(0)\A b ( Pl ,s)) 



(3) 



The amplitude T^ v is decomposed into 14 invariant form-factors 7$, Si analogous to the 
Wi,Gi in (0) and can be computed by using HQET methods p fil^fi^^JTOl , ^ . Keeping 
terms up to 1 jm\ one obtains 



1 2m,, 

(m b - v .q)(l + X b ) + —±(K b + G b 



(v.qf 



2A 



2An 



A 2 
^0 



where 



m b ( K b + G b ) _ m 2 b G 
2A 3A 5 



+ ^ ^ttK-m). 



T2 = ^ + X b ) + —(K b + G b ) + 

m b (K b + G b ) Am 2 b K b v.q 2mjjG b 
A 3A§ + 3A 2 ' 



2 ' 



^4 = ^2(^6 + ^&)> 

-1 , s 2m b 2m^^ m|Gfe 

Ts = 2A^ (1 + Xb) ~ — {Kb + Gfc)_ Ag - + l4~ 
1 + e 6 5m 6 km\K b( 2 2 

^-^aF' 



3 

^4 = 0, 

-2m b K b 

b5 ~ 3A 2 ' 

m fe (l + e fe ) 5m b if fe 5m 2 b ^m\K b 2 2 

<J6 = : : r^v.qK b H 7-5— (0 — (f .0) ) , 

2A 6A 3A 2 , y 3Ag v v ' ' 

l + e b (2m b + 3v.q)m b K b Am 2 b K b 2 2 

57 = -2AV + 3Al 3AjT 

5 8 = mbi \ +eb) + ^ + ^t.^ 6 - ^(g 2 - (v.qf) , 

2A 6A 3A 2 , H 3Ajj w v y; ; ' 

2A 3A 2 , 3Ag w v H! ' ' v 7 



-2(m 6 - v.q)m b (K b + G b ) 8m b K b 2 2 2m^ 6 , . 

x&= a 3AT (g + (5) 

We use the notation of |TIJ throughout the paper. K b is related to the mean kinetic energy 
of the heavy quark in the A b baryon 

(iD) 2 u 2 
K b = -Y t (Ab(p,s)\b v (x)^-Lb v {x)\A b (jp,s)) = « 0.013, (6) 

s Zm b Zm b 

where we used fi 2 ~ 0.6 GeV 2 |18j and m b = 4.8 GeV. The spin-dependent contribution 
e b is defined by 
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(A b (p,s)\h X J 5 b\A b (p,s)) = (l + e b )s x . (7) 



with e b ~ — \K b . An estimate of the spin-dependent parameter e b has been given in [IK 



with the result e b = —^K b , based on an assumption that the contribution of terms arising 
from double insertions of the chromomagnetic operator can be neglected. A zero recoil 
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sum rule analysis gives the constraint e b < —%K b [0 which puts the estimate of [I9~| at 



the upper boundary of the constraint. We use the value of [19] keeping in mind that the 
numerical value of e b could be reduced in more realistic calculations. 

In order to make our presentation as complete as possible we have retained the 
chromomagnetic contribution proportional to G b in Eq.@ although it is zero for the A b 
system, i.e. 

G b := Y,(MP,s)\b v (x) ( ~ 9 ^f P ) b v (x)\A b (p,s)) = . (8) 

The reason is that the general helicity formalism introduced later on can also be applied 
to B meson and Q b baryon decays where G b ^ 0. 

The denominator factor Ao is given by 

A = (m b v — q) 2 — m 2 + ie (9) 

The imaginary parts of inverse powers of A , which are needed for obtaining the structure 
functions Wi,Gi, can be obtained with the help of 



Im 
Im 





= ~ n 6 
2m b 


Ajj) 


_7F 5 
Ami 




— 71 




16m| 



-m 2 + q 2 + ml \ \ 
2m~ b J J 
-m 2 + q 2 + m 2 b Y 
2m b ) i 
-m 2 + q 2 + ml * 



Various subsets of @ have appeared in the literature before ]]9| JT0| , [T6|1 . We have recalcu- 
lated them and collected them together for ease of reference. 
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III. HELICITY STRUCTURE FUNCTIONS AND ANGULAR DECAY 

DISTRIBUTIONS 



The hadronic structure for the transitions A&(s) — > Xj is fully specified by the 
absorptive parts of the 14 structure functions listed in Eq. (|4]). In order to obtain the 
full decay distribution for the inclusive decay A b {s) Xj + l~(si) + V\ one needs to 
contract the hadronic tensor W^ v with the known leptonic tensor L^ v . Traditionally the 
contraction L^W^ is done in covariant fashion. Here we advocate a different approach 
and use helicity techniques to write down the relevant decay distributions. The advantage 
is that the angular decay distributions involving helicity angles are given by simple linear 
combinations of the helicity structure functions. The use of helicity techniques to describe 
angular decay distributions in exclusive semileptonic decays is widespread by now 
and is easily generalized to inclusive semileptonic decays. 

Let us begin by writing down the relation between the full spin-dependent set of 
14 helicity structure functions and the set of 14 invariant structure functions in Eq. (|2]). 
One has 

W++ = {ql - q 2 )G 3 - (d + W 3 )p - q G 9 -G 8 + W x , 

W++ = ~{{ql - q 2 )G 3 - {G 1 - W 3 )p - q G 9 - G 8 - W x ) , 

W++ = ~{{q 2 - q 2 )G 3 + (G x - W 3 )p - q G 9 - G 8 - W x ) , 

WZZ = {q 2 - q 2 )G 3 + (G 1 + W 3 )p - q G 9 - G 8 + W x , 

{-{{2{q 2 G 5 - G 6 )q - {G x + 2G 7 )q 2 + q 2 G 2 + q A G 4 )p - q 2 W 2 - 2q q 2 W 5 - q 4 W A + q 2 W x )) 



tt — 2 



{{2{q 2 G 5 - G 6 )q - {Gx + 2G 7 )q 2 + q 2 G 2 + q 4 G A )p + q 2 W 2 + 2q q 2 W 5 + q A W A - q 2 W x ) 



tt ~ q2 



W, 



{ql ~ q 2 )W 2 + {2q G 6 - g 2 Gi)p - p 3 G 2 + q 2 W x 



00 — 2 

q l 



{ql ~ q 2 )W 2 - {2q G 6 - q 2 G x )p + p 5 G 2 + q 2 W x 



oo " g2 



++ 



pq W 2 + pq 2 W 5 - q$G 2 - q 2 q 2 G 5 + 2q 2 G & + q q 2 G 2 + q q 2 G 7 + q 4 G 5 - q 2 G e 



ot ~ q2 



W { 



pq W 2 + pq 2 W 5 + qlG 2 + q 2 q 2 G 5 - 2g 2 G 6 - q q 2 G 2 - q q 2 G 7 - q A G 5 + q 2 G 6 



ot ~ q2 



-2{pG 6 + g G s + q 2 G 9 
V¥V2 
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w+ _ _ -2(pG 8 + q G 6 + q 2 G 1 





V¥V2 


-2(pG s 


- q G 6 - q 2 G 7 ) 






2(pG 6 - 


q G 8 - q 2 G 9 ) 



W3T= " ° fey ' ■ ("J 

V? 2 V2 



Aa,A' a 



We denoted here p = \Jq 2 — q 2 - The helicity structure functions W AwA/ 6 are defined by 
= (2vr) 3 ^5 4 ( Pl -g-^)(X J |Jf|A b ,A A J e ;(AH/)(A ; , ) A;jj; t |X,) ei ,(AV). (12) 



Here Aa 6 = ±1/2 is the helicity of the A b and Avk = 0, ±1, t are the helicities of the virtual 
J¥-boson (spatial: Xw = 0, ±1; temporal: \w = t). Note that there are no zeroth order 
parton model and kinetic energy contributions to the structure functions Wtt and WZZ 
because of angular momentum conservation. We shall return to this point later in this 
section when we discuss the contribution of the spin-dependent matrix element e b . Note 
that from Eq. (|T2"|) one has the hermiticity relation 

Aa, A', * A'. Aa, 

W X \, A » =Wy\ b . (13) 

Since the helicity structure functions are real in our case one can drop the complex 
conjugation sign in Eq. (p~3|) . The helicity structure functions are defined in the A b rest 
system. We therefore need to specify a z-axis which we take to be along p x (see Fig.l). 

As noted before the full angular decay distribution of A b (s) — > Xj + l~(si) +i>i includ- 
ing all polarization effects is completely determined by the set of 14 helicity structure func- 
tions. The necessary manipulations involving Wigner's D J mm , functions are standard and 
well documented in the literature [pl] , p2]| (see also Sec. IV). Here we closely follow the pre- 
sentation of @,§,0]. For example, for the five-fold decay distribution in q , q 2 , cos 6, cos Qp 
and (p into negative (dr - ) and positive (dr + ) helicity leptons we obtain ' 



dr" 2G 2 \V b] \\q 2 -m 2 ) 2 Jq 2 -q 



dg dg 2 dcos Gdcos 9pd0 3(27r) 4 g 2 



^Similar decay distributions have been written down in [23], where the 0(a s ) corrections to 
unpolarized t — ► b decays were evaluated. 
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P++ iwii + w: 



[wzz + w- 



+7 (*+ K+ + 



++ 



+ p- (w~- 



w 



p+- [WI + Wfc J sin 9 cos 



w 



;i + cos 2 e) 



cos 9 



2^ 

+ 7^|P+- (W+f " sin26cos0 

dr+ 2G 2 |H,f(g 2 -m2) 2 v /g 2 -g2 



(14) 



dgodg 2 d cos Od cos 9 p( 



3(2vr)V 



2g 2 



( P++ (w±t + w++) + P __ (wrr + w--)) - sin 2 e 

3 



P++WJo + + p—Woo ) -cos^G 



+ 3 (p ++ ^ t + + p— WoTj cos 

- ^v/2p+- (W+f - W+-) sin 9 cos0 

- -V2p+- (W+ - sin 26 cos< 



(15) 



For A c — > X s + l + + z/; decays one has to effect the replacement dr T <-» dr ± and one has to 
change the signs of the parity violating contributions proportional to cos and sin sin <j) 
in Eq.([TJ]). The differential rate into unpolarized leptons is simply dr + +dr~. The polar 
angles and 0p and the azimuthal angle are defined in Fig.[l]. We have rotated the 
density matrix of the A& to the z-axis such that one has 

, / 



p Ab (cos0 P ) 



l + Pcos0 P Psin0 P 
Psin0 P l-Pcos0 P 



(16) 



Note that for A^'s from Z-decays one expects that the A&'s are longitudinally polarized 
with backward polarization. Thus, for A^s from Z-decays, the direction of P in Fig.|I] 
coincides with the boost direction that brings A& to rest (P > 0). 

The longitudinal polarization of the r-lepton is given by 

dT + - dT~ 



PL 



dr+ + dr- 



(17) 



9 



We emphasize that the longitudinal polarization of the r calculated from flT4] , |15| ) and 
(|17|) refers to the (r, v T ) cm. frame which differs from the longitudinal polarization of the 
t in the A& rest system as calculated e.g. in |||T7). A slight adaptation of the lepton-side 
density matrix elements in flU, |T5|) as described in M will yield the longitudinal polar- 



ization in the A& rest system. Put in a different language the two respective polarizations 
are related to one another by a Wigner rotation (see e.g. ||24|). A direct computation 
of the longitudinal polarization of the r in the A& rest frame including correlation effects 
will be presented in Sec.V. 

The quasi three-body decay A&(s) — > Xj + l~ + i>i is described by three kinematic 
invariants. A particularly convenient choice is the one used in Eqs.flTJJ, ITS]) in terms of 
q , q 2 and cos 6. Working in (g Q , q 2 , cos O) phase space has big technical advantages as can 
be seen in the following. The integration over cos is trivial. The integration over g is 
almost trivial as can be seen by the following reasoning. The leading order parton model 
or free quark decay contributions are proportional to the 5-function and the go-integration 
amounts to the substitution q = b 2m J — in these terms. The corrections to the parton 
model contributions involve derivatives of the 5-function. Using partial integration the 
derivatives can easily be shifted to the integrand functions without encountering surface 
term contributions because, to the requisite order in l/m^, the two-dimensional parton 
phase space never touches the boundary of the three-dimensional particle phase space, 
except at maximal q 2 and go where the integration measure is zero. The only nontrivial 
phase-space integration that remains to be done is with regard to q 2 . 

If desired, the transformation to the usual (go, g 2 , Ei) set of variables can be done 
with the help of the relation 

g (g 2 + m 2 ) -2q 2 Ei 



cos 6 



g 2 - q 2 {q 2 - m\ 



However, as can be easily appreciated by substituting Eqs.(|T^) in the differential rate 



Eq.(14, |T5|), the integration of the corresponding differential rate function is more cum- 
bersome for the following two reasons. First, when doing the go-integration, one has to 
carefully consider the contributions from the surface terms induced by the derivatives 
of the 5-function from the OPE expansion (see Eqs.(|K])). These may lead to the ap- 
pearance of spurious singularities when calculating polarization type observables Plffl. 
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Second, one remains with two nontrivial (Ei and g 2 ) integrations, apart from having to 
carefully consider surface term contributions when doing the go-integration. 





6 — > c 


b — > u 


2 / 2 

77 = m T /mi 


7] = 


1 1 2 

rj = m T J mi 


T] = 




0.0389 
(0.0376) 


0.177 
(0.172) 


0.133 
(0.124) 


0.352 
(0.333) 


fz 


0.0330 
(0.0352) 


0.332 
(0.344) 


0.124 
(0.136) 


0.635 
(0.666) 


fp 


-0.0185 
(-0.0208) 


-0.105 
(-0.112) 


-0.093 
(-0.124) 


-0.282 
(-0.333) 


fr 


0.0183 
(0.0208) 


0.104 
(0.112) 


0.092 
(0.124) 


0.279 
(0.333) 




-0.0282 
(-0.0306) 


-0.322 
(-0.331) 


-0.116 
(-0.136) 


-0.629 
(-0.666) 




-0.0385 
(-0.0376) 


-0.176 
(-0.172) 


-0.132 
(-0.124) 


-0.349 
(-0.333) 


ff- 


-0.0397 
(-0.0349) 


-0.208 
(-0.213) 


-0.123 
(-0.128) 


-0.418 
(-0.431) 




-0.0239 
(-0.0263) 


-0.179 
(-0.186) 


-0.103 
(-0.128) 


-0.390 
(-0.431) 



Table 1. Numerical values for polarized and unpolarized helicity structure 
functions into negative helicity leptons. Column 2: b — > c, I — r; column 
3: b —>■ c,l = e (m e = 0); column 4: b — » u (m u — 0), I — r; column 5: 
b —>■ u,l = e (m e = 0). 



Returning to the differential rate Eq.(14, [15|) one first does the go-integration which, 
after integration by parts, amounts to a mere substitution as noted before. The results 
of the go-integration are listed in Appendix A. Next one integrates over q 2 in the limits 
mj < q 2 < (mi — m 2 ) 2 . Numerical and analytical results of the g 2 -integration can be 
found in Tables 1 and 2 and in Appendix B. 
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b — > u 




0.0081 
(0.0079) 


0.0198 
(0.0186) 




0.0076 
(0.0080) 


0.0218 
(0.0234) 




0.0079 
(0.0080) 


0.0231 
(0.0234) 




0.0066 
(0.0072) 


0.0210 
(0.0234) 




0.0041 
(0.0046) 


0.0154 
(0.0186) 




-0.0067 
(-0.0072) 


-0.0213 
(-0.0234) 




-0.0066 
(-0.0072) 


-0.0208 
(-0.0234) 




-0.0077 
(-0.0080) 


-0.0222 
(-0.0234) 


1 ST 


-0.0077 
(-0.0076) 


-0.0209 
(-0.0205) 




-0.0055 
(-0.0060) 


-0.0181 
(-0.0205) 



Table 2. Numerical values for polarized and unpolarized helicity structure 
functions into positive helicity r-leptons for b — > c and b — > u transitions. 
Parameter values as in Table 1. 

We use in the Tables 1,2 the following numerical values: m T = 1.777 GeV, m c = 
1.451 GeV, m b = 4.808 GeV. K b is given by K b = f4/(2m 2 b ) with /4 = 0.6 GeV 2 and for 
e b we adopt the value e b = —2K b /3 as explained above. The values in brackets show the 
free quark decay (FQD) results. 

In order to simplify our notation we introduce a set of unpolarized and polarized 
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reduced differential rate functions dt~ and dtf~, respectively, for the decays into negative 
helicity leptons. We define scaled variables q 2 = q 2 /ml, % = qo/rrib, P = rn 2 /ml and 
rj = rn^/m 2 and write 

= iQ^-^mtt + wzz + {-){w-- + wzz)) 

^ = 16^£/(^ + o + + (-)WS-) 

rlf (P) ~ (a 2 - n) 2 

dip = 16 W iyj J(wg _ + _ 
dg^ g 2 

M: = i 6 (i!_# /(in - _ ( i9) 

where df stands for either df y or df p ~ etc. with the corresponding signs specified 
on the r.h.s. of Eq. (|T9|). Accordingly we define reduced differential rate functions for 
decays into positive helicity leptons. 

( p )+ „ aH p )~ 



u 



dq 2 2q 2 dq 2 



dff )+ n dt[ p) - 



dq 2 2q 2 dq 2 
dff )+ rj (q 2 -v) 2 



dg z 2q z q z 

+ HWi 



df£ )+ Of-?/) 2 



df p+ _ n df p ~ 



(20) 



dg 2 2q 2 dq 2 

Note that the structure function combination I{WZf +Wf+~) does not appear in Eqs. ([19|) 
and (|20|) . This combination can only be measured through the transverse spin components 
of the r-lepton as discussed in the next section. The differential rate into negative helicity 
leptons can then be written as 

dr- T b 



dg 2 dcos Odcos 6pd0 An 



x 
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where 



3 (dTTj n ^ dH- 



3 dT] 



+- 



+ 



4 \ dq 2 
3 fdTz 



+ PcosGp 



4 \ dq 2 
3 rfff 



+ Pcos6 P 



dq 2 
dH 



;i + co S 2 e) 



sin 2 e 



cos 



2^/2 dq 2 



-P sin 6p sin 6 cos — 



3 dY p A 
4^2 dq 2 



-P sin 6 p sin 26 cos (j> 



(21) 



\V b3 \ 2 G 2 F ml 

1927T 3 



(22) 



For the rate into positive helicity leptons we have 



= — x 



dg 2 dcos6dcos0pd0 Air 
3 



4 V dq 



u + PcosGp- 



dq 



fa 



sin 2 



+ 



2 V dq 



L + PcosQp^r- I cos 2 



^2 



<ig 



£2 



3 
+ 2 



+3 



L + PcosGp- * 



<ig 



$2 



<ig 



^2 



cos 



+ - V2—^P sin 6p sin 6 cos + - V2—4—P sin 9 P sin 20 cos 



dq 



dq 



(23) 



For quick reference we shall adopt a generic labelling for the various helicity structure 
functions, namely we write for dr[P~/dg 2 , etc. An inspection of the g 2 dependence 

of the helicity structure functions shows that the structure functions , , F^~, 
f/( p ) + , , S'( p ) + and SL^ pS)+ can be integrated analytically. In contrast to this, the 

integrations of the structure functions A^" , ST^ pS)+ and A^ p ) + lead to incomplete 

elliptic functions and will therefore be performed numerically. 

In Appendix B we list analytic results for the totally integrated unpolarized structure 
functions U~, L~ , F~, U + , L + and S + as well as for the total rate function proportional 
to (U~ + U + + L~ + L + + 3S + ). Numerical results are listed in Tables 1 and 2, where we 
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have used the following set of input values: m b = 4.8 GeV, m c = 1.45 GeV, m u = GeV, 
m T = 1.777 GeV, K b = /4/(2mg), /4 = 0.6 GeV 2 and e b = -2K b /3. In order to assess 
the importance of the nonperturbative contributions we also list the zeroth order parton 
model values (K b = e b = 0) in brackets. Judging from the numerical entries in Tables 1 
and 2 the effect of the nonperturbative contributions can go both ways, i.e. depending on 
the particular structure function, they can enhance or decrease the parton model results. 
The nonperturbative contributions are generally small, at the few percent level. Notable 
are the structure functions U p ~ , A p ~ , SL + , U p+ and L p+ where the nonperturbative 
contributions exceed 10%. 

Fully integrated values of the longitudinal polarization of the r can easily be con- 



structed in analogy to Eq. (17). It is also evident that there are no transitions into 
positive helicities for mi = which explains why Table 2 has only two columns. 

The structure functions combinations (U~ + F p ~ ± (U p ~ + F~)) are of particular 
interest since they are directly proportional to the nonperturbative spin parameter e b . 
They receive contributions from the helicity ±3/2 configurations, which are neither pop- 
ulated by the parton model contribution nor by the spin neutral kinetic energy term K b . 
In fact, for the rate combinations (U~ + F p ~) and (U p ~ + F~) one finds 

f^ + f^ = -e b frj(^ = 0) 

m 

+8 



1 - 7t] - 7r] 2 + rf - 7p + YlT) P - 7t] 2 P - 7p 2 - 7t] P 2 + p A ) (24) 
p 2 (r] -l)ln + 7] 2 {p 2 - l)ln ' 



2^ / V 2 ^ 



and 

T p ~ + T- = -e b Tp(K b = 0) 
= -e b 



(-1 + ?7 + 2VP - p)(l — 77? — 7 V 2 + v A + 2^p- 12 V ^p - 2 V 2 ^p - 17 p 
+38r/p - 7ifp + 28v^ 3 - 1277VP 3 - 17p 2 - 7i 1 p 2 + 2^/f + p 3 )/3 

(25) 



+4r7 2 (l-p) 2 ln' 7 



where the analytical results on the r.h.s. of Eqs. ( p^j) and (|25|) can be obtained from the 
closed form expressions for and Tp in Appendix B by setting K b = 0. 
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Numerically one has 

f £ + tfr = -0.0376 • e b 

f + t F = 0.0208 • e b (26) 

It is clear that the contribution of K b could be extracted from taking appropriate 
linear combinations of the fully integrated structure functions listed in Appendix B. Since 
the coefficients needed in this extraction involve high powers of the masses in the process, 
which are uncertain, it would be desirable to find an observable directly proportional to 
K b . Such an observable appears in the measurement of the transverse polarization of the 
r as discussed in the next section. 



IV. TRANSVERSE POLARIZATION OF THE r-LEPTON 

In this section we present results on the transverse polarization components of the 
r-lepton in the (r, i> T ) cm. frame. The two transverse components are conventionally 
divided into the transverse perpendicular component in the lepton-hadron plane and the 
transverse normal component out of the lepton-hadron plane. The latter component is 
not affected by the boost from the A b rest frame to the (r, v T ) cm. frame and, after the 
appropriate integrations, can therefore be compared with the corresponding calculation 
done in the A b rest frame in ||. 

The 2x2 density matrix of the r-lepton can be obtained from the master formula 



V — {^(ei,e P )} Alii , (27) 



dg dg 2 dcos 6dcos 9pd0 3(2ir) 4 q 2 
where @ 

{W(6,0,0p)} A;A; = A £ [(-l) J+J 'h Xl ihlne-^-™ ') 

m,m' ,J,J' ,X b ,\' b 



(2f 



The helicity amplitudes h Xi ± (A; = ±|) for the decay W fj_ shell — > r + v T appearing in 
the master formula are given by || 
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1 1 

' 2 2 



hii 

2 2 



8(g 2 - m 2 T ) = 2V2r 



i rat. , 
— ~h 1 1 
2q 2 



-h_ii . 
2q 2 2 2 



(29) 



The diagonal elements of the density matrix {M^}a ; aj have already been written down 
in the main text. Here we list the nondiagonal density matrix element relevant for the 
transverse polarization components of the r. 

For the unnormalized transverse perpendicular component of the polarization vector 
one needs to calculate {W} x = {W} + ^ + One has 



{w}' = {w{e,t/>,e P )} 

3 



{W(Q,<t>,Q P )}- 



2s/2 
3 

~4\71 



{ P++ mt-w++)+ P -{wi 



M -mi) 



( P++ (w^ + wl+ 



2W++) 



+ p__ (W~- + WZZ - 2 W o - )) sin20 



-—(p ++ W+ + + p __W - t 



-2P+-( W <H 



+- 



w: 



sin G 

cos 26 cos cb 



W+n + W+- - Wir) cos 9 cos . 



w+ t -)cos< 



(30) 



Note that the 14 th structure function combination {W^T + W^f) that was missing from 
the rate expressions in the main text makes its first appearance in the transverse normal 
polarization. As Eq.(^) shows, A^-polarization as well as a determination of the r- 
lepton's transverse polarization is necessary for a determination of the complete set of 14 
structure functions. 

The unnormalized transverse normal component is given by the combination 
i(W+- — W- + ). One obtains 



{wy = i({w(e, 0, e P )} + _ - {w(e, <f>, o P )}_ 



m 



= (q 2 -m 2 T ) 



w: 



w t + + ~ 



w: 



cos O sin 



t+ ~ W -t ) sin ' 



(31) 
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Note again the contribution from the 14 th structure function combination (W 7 ^ - + W^f ). 
It is quite remarkable that the two combinations of helicity structure functions appearing 
in the transverse normal polarization of the r in Eq. fl3l|) can be seen to be entirely 
determined by the nonperturbative contribution. We have no simple explanation of 
this fact. 

The go-integration of the relevant density matrix elements fl2"T|) can easily be done 
as described in the main text. In fact, the relevant go-integration of the 14* h structure 
function combination is listed in Appendix A. The remaining g 2 -integration is then done 
numerically. The numerical result is presented for the case b — > c using the same numerical 
values as in Sec. III. We obtain 

dT x r 

— — — - = — [-0.03777 sin 6 + 0.00725 sin 26 

dcos Bdcos Bpd<p An 

+P cos P (O.O1519 sin + 0.01907 sin 26) 

+P sin Bp cos 0(0.01682 + 0.00091 cos 6 + 0.01714 cos 20)] (32) 
U b (-0.000908 + 0.000323 cos Q)P sin Q P sin (33) 



dcos Gdcos 9pd0 An 



The transverse perpendicular polarization is dominated by the zeroth order parton con- 
tribution with angular coefficients comparable to the entries in Tables 1 and 2. The 
transverse normal angular coefficients are quite small as expected since they are propor- 
tional to Kb. 

As mentioned before the transverse normal polarization is not affected by the boost 
from the lab frame to the (t~v t ) cm. frame. In order to compare our results with 
the corresponding results in Ref.|| we integrate the transverse normal polarization with 
respect to cos0 and cos Op to obtain 

dlw = -£ i2rf ^/5 /(H/ - + w -°- + w >+ - < 34) 

As noted above, the linear combination I(W^T + W^q + W^T — W^f) is proportional to 
Kb only and has no zeroth order partonic contribution. The g 2 -integration is easily done 
and one finds 

^ = —PAsm<f) (35) 
d0 2tt y y ' 
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where 



A = 2nK b ^q 

-4^(1 - 2p + rjp + p 2 ) In 
— 4p(l — 2i] + r/p + rf) In 



V^R(-2 - 5r/ + rf - 5p + p 2 + 10r/p)/3 
1 +r/ - p - y/R\ 



2 VV J 
A -y + p- \fR" 

2VP 



(36) 



This result agrees with the result in II when their corresponding y-distribution is inte- 
grated with respect to y. Numerically one has 



A = -O.lLK* 



(37) 



The results on the r polarization discussed in the last two sections refer to the 
(r, v T ) rest frame which, because of the elusiveness of the neutrino z/ r , may not always be 
easy to construct. For some applications it may be preferable to avail of the longitudinal 
polarization of the r in the A^, rest frame. This is the subject of the next section. 



V. POLARIZATION CORRELATIONS IN THE A B REST FRAME 

In this section we determine the polarization of the r in the At, rest frame using a 
calculational technique originally proposed in |26| which involves an averaging over the 
azimuthal angle <fr (see Fig.l). In this way we determine the polarization components of 
the r in the plane spanned by the r and the polarization vector of the A&. 

The method of |26| is based on applying the unitarity relation to the amplitude for 
forward fA& — > fAb scattering (see Fig. 2) 

T(s) = i(A b {v, s)f{p T )\ J dxTHlyixWwWlMv, s)f(p T )) , (38) 

where 

H w = 2V2V jb G F \j^P L b][e^P L u e ] (P L = i(l - 75 )) (39) 

is the interaction Hamiltonian responsible for the semileptonic decay b — > jr^u T and 
s = (m^ b v + p T ) 2 . By inserting a complete set of states in ( |3"5| ) one obtains 
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JmT(s) = l J2 Jdfi(X J )dfi(iy T )(2n) 4 5 4 (m Ab v + p T -p Xj -p 1/r ) (40) 

Xj,v r 

x |(X^ T |^(0)|A 6 (i;, S )f(p r ))| 2 . 



The phase-space volume element is d/i = ^yi 2E ■ 

Comparing (|40|) with the inclusive semileptonic decay rate 



dr = \{X J v T {p UT )r{p T )\n w \X b {v,s))\ 2 {2-K f5\m Ab v - p T - p Xj - ^Jd/i(X J )d yU (z/ T )d yU (r) (41) 

one obtains the final formula relating fl38|) to quantities of experimental interest: 

dr = -i-l^ld^d^Im T(s) . (42) 

Our problem is thus reduced to computing the function Im T(s). This can be 
done with the help of an operator-product expansion combined with the heavy mass 
expansion as discussed in Sec. II. First the heavy hadron is replaced with a heavy quark 
with momentum pt, = m^v + k. From the 1-loop diagram in Fig. 3 one reads off the 
lowest-order expression for T(s) 

T(s) = 8G 2 F \V 3b \ 2 [u{p T ) lfll a lv P L u(p T )} [u(o)7yfPL«W] Iaflik) (43) 
where we have defined 

a/3 J (27r) rt \(p T — q) 2 + ie][(mbV + k — q) 2 — m 2 + it] 

We will be interested in polarized r leptons in the final state. Thus the leptonic 
amplitude in (|43| ) should be replaced by 

uiPrh^VMPr) - Tr {^(A + + 75^)7m7 q 7^l} , (45) 

with s T the spin vector of the r lepton. 

The integral I a p(k) in ([44]) can be calculated by combining the denominators with 
a Feynman parameter x. The result is a complex function with an imaginary part given 
by 
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-Urn I a/3 (k) = (46) 

Z7T 

where 

s = x |— m^(l — x) + — (m(,t> + /c) 2 (l — x) + 2(1 — s)p r ■ (m^f + . (47) 

In (|46|) we have denoted xi the root of the equation s(x\) = 0. It is given by 

x x = \ L _ (48) 

1+77-2/1 + 

= x + — I 2k ■ (v p T ) + k — H , 

1+rj-y \ l + r}-y J 

where x = 1 — p/(l + ry — y) is the value of Xi for = 0. We have introduced reduced 
momenta k = k/m b , p T = p T /m b and have expanded in powers of k up to second order. 

The integration in fl4"6| ) can easily be performed with the result 

9aj3 



Im I aP {k) = ^\\px\ + I (y - y) H " ^ + ^ + 2pT -i v + ~ k ) 



- ^ (j ~ + k - PrUv + fc - Prh , (49) 

where X\ has to be replaced with the expanded form of Eq. (|48|) . 

In physical applications we are interested in heavy hadron decay, rather than free 
heavy quark decay. One can obtain the corresponding scattering amplitude T(s) by 
replacing the hadronic spinor expression in ( f4~3|) with expectation values of the appropriate 
operators 

u(v)f(k a )u(v) - (A b (v, s)\bf (^) b\A b (v, s)) . (50) 

These matrix elements can in turn be expanded in powers of l/m b with the help of heavy 
quark effective theory (HQET) methods as discussed earlier. Referring to |i"3| , |l0f for 
calculational details, we only give the final substitution rules needed for computing Im 
T(s). 
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a) terms of order k 



b) terms of order k l 



7 M -> *V (51) 
7 M 75 -> -v(l + e 6 ) (52) 



fc M 7^ -> ^fe(-2^ + 5w M ^) (53) 

2 

V^75 K b (Vf,s u + -v u s^) (54) 

c) terms of order A; 2 

~ ~ 2 

&/A7a -g-MfiV - v^v u )v a (55) 
2 

fc M ^7a75 ~> -^ K b{g^u ~ Vfj,V v )S a (56) 

It is now a simple matter to combine (|3S|J3I3) and (^T]-^) into ( f4"2"D and extract the 
decay rate. It can conveniently be written by splitting the rate into two terms 

dr = dr + dP> = 1 (dT(s T ) + dT(-s T )) + £ (dT(sr) - dT(-s T )) • (57) 

dr and dr p represent the decay rates into unpolarized and polarized leptons. For the 
unpolarized rate we obtain 

1 dr 



T b dyd cos 6 T 

where 



^-4 V (A(y) + (p T -s)B(y)) , (5* 



A(y) = x 2 (-3y 2 + 6y(l + rj) - 12 V ) + x 3 (y 2 - 3y(l + rj) + Srj) (59) 

+ K b {2[(1 + r?) 2 - y 2 } - Ax [y 2 - y(l + rj) + 2(1 + r, 2 )} 

4 32 
+ a; 2 [4|/ 2 - %(1 + rj) + I677 + 10(1 + r/ 2 )] + x 3 [--y 2 + 4y{l + rj)-—r)- 4(1 + r, 2 )} 

_ ^l-yni + vKl-Xoni-Sxp) + 2 (l-#(l-a:o) 2 (l-4x ) | 
l + r]-y (l + r]-y) 2 J 

= {(1 + e b )[6x 2 (y - 2rj) - 2x 3 (l + y - 3rj) 



A(y + 2) - 8x (2 + r, - y) + 8x 2 (l + 2r/ - y) + -a^(y - 877) 



i (l- V )(l-x ) 2 [3 + V -2x ( V + 2)} 4 (1 - 77) 3 (1 - x ) 2 (l - 4x ) 



l + 77-y {1 + 7} -yf 
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We have defined 6 T as the angle between the A b spin and the r lepton momentum direction. 
For the decay rate into polarized leptons dr p we obtain 

1 dr? 



Tb dyd cos 9 r 
where 



^Jy 2 -4vVv{(v ■ Sr)(p T ■ s)A*{y) + (v ■ s T )B p (y) + (s ■ s T )C p (y)} , (60) 



A"(2/) = -24(l + e 6 )l| 



■ I n 



-4K h 



2 + 2x - 5xq + -a?o 



(61) 



2 (l+r ? )(l-x ) 2 (2-x ) 2 (1- r?) 2 (l - x ) 2 (l - Ax ) 



-y 



[1 + 7] -y) 2 



B p (y) = -6x 2 (2 -y) + 2x^(3 - y - rf) - 4K b [-y - 2 V + 2x (l + 2 V - y) (62) 

+ 2xl( - 2 -, + y) + * 3(3 - ,) + (l-,)(l-^[-l-3„ + 2. o(1 + 2„)] 

3 1 + 7] — y 



+ 



(l-77) 3 (l-a;o) 2 (l-4xo) 

(i +v - y) 2 



C p {y) = 2x 3 (l + e b )(l + r]-y)- 2K b 2x (l + tj + y) + x 2 (-4 -4>q + y) 



(63) 



+xl(2 + 2 V - -y) 



2x (l-i 1 ) 2 (l-x Q f 
l + r]-y 



We have checked that this formula agrees with the various particular cases pre- 
sented in the literature. Thus, |17| compute the longitudinal polarization asymmetry 
corresponding to an unpolarized decaying baryon and [[| give the full results for vari- 
ous lepton polarizations from the decay of polarized A& baryons. The method presented 
here has the advantage of performing the integration over the neutrino phase space au- 
tomatically. Unfortunately, this aspect can prove to be also a limitation: because of 
integrating over all possible neutrino momenta, all information about the decay plane is 
lost. Therefore this method can only be applied for obtaining lepton polarization asym- 
metries averaged over the position of the decay plane, i.e. after azimuthal averaging. The 
azimuthal averaging does not affect the longitudinal polarization component and thus our 
results can directly be compared to corresponding results in the literature [§,[L7j] . We have 
also checked that we agree with on the transverse polarization component in the plane 
spanned by the p T and P which can be obtained from with the appropriate azimuthal 
averaging. The transverse polarization component normal to this plane averages out when 
doing the azimuthal averaging. Our results for the transverse normal component of the 
r-polarization have been presented in Sec. IV. 



23 



VI. SUMMARY AND CONCLUSIONS 



We have analyzed the inclusive semileptonic decays of polarized A b baryons into 
polarized r-leptons. We discussed spin-spin, spin-momentum and momentum-momentum 
correlations between the spins of the A& and the r, and the momenta of the virtual W 
(or recoil momenta px) and the r. Using helicity techniques we presented detailed results 
on the above angular correlations involving a three-fold angular decay distribution in the 
three helicity angles that can be defined for the process. By taking suitable combinations 
of helicity structure functions we identified observables that are directly proportional to 
the contributions of the 0{l/ml) nonperturbative matrix elements. In the helicity method 
one determines the r-polarization in the (t~, v T ) rest frame. We give also results on the 
r-polarization in the A^-rest frame using an elegant loop calculation. 
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APPENDIX A: Qo-INTEGRATED STRUCTURE FUNCTIONS 

In this appendix we list the results of integrating the structure functions in Eq.flllD 
with respect to go- Our results are given in terms of the integrals 



where W stands for any of the linear combinations of the helicity structure functions in 
Eq.(|ID. One has 

I(W++ + + WZZ + WZZ) = (1 - K b )p{-q 2 + p + l)/2 + AK b p/3 

+G 6 (p 2 (15(g 2 - p) - 11) + g 2 (3g 2 - 3p) - 7) - 4p + 4))/(6p) 
I(W++ + WZZ - WZZ - WZZ) = (1 + e b )p 2 + 2K b {-2p 2 - q 2 )/3) 

I(W+ + + W 00 -) = (1 - K b )p(Ap 2 - g 4 + q 2 p + g 2 )/(4g 2 ) - AK b p/3 
+G b (p 2 (15(-4p 2 + g 4 - q 2 p) - 59<f - Yip + 12) 
+ <f (<f (3(<f - p) - 7) - 4p + 4)) /(12pg 2 ) 
J(W ++ - WV -) = (1 + - l)/(2g 2 ) + K b (-2p 2 - p + l)/3 

- + - Wrr) = -p 2 + 2K 6 (2p 2 + g 2 )/3 + G b (Wp 2 + 2q 2 + 3p - 3)/3 

I(W++ - - WZZ + WCl) = (1 + e 6 )p(g 2 - p - l)/2 + K b p{-3q 2 + 3p- 5)/6 

/(W+r + Wo+) = (1 + e b )p(g 2 + p - 1)/(2V2^) 

+K b p(q 2 -p+l)/(3V2^f) 
I(W+Z - W+-) = -(1 + e b )p 2 /(v / 2 v^) + ir 6 (-p 2 + 2q 2 )/(3V2^P) 
KW++ + ^--) = (i _ K b )p(4p 2 - g 4 + g 2 p + g 2 )/(4g 2 ) 

(p 2 (5(-4p 2 + g 4 - 5g 2 p) - 9g 2 - 4p + 4) 
+ g 4 (g 2 -p-l))/(4pg 2 ) 
I(W+ + - W t r) = (1 + e fc )p 2 (p - l)/(2g 2 ) 

+K b (-2p 2 p + 2p 2 - g 2 p + g 2 )/(3g 2 ) 
I{W ++ + W -) = f(_ p + i)/( 2 f ) + ^( 2 p 2 p - 2p 2 + g 2 p - g 2 )/(3g 2 ) 
+G 6 (10p 2 p - 2p 2 + 5g 2 p - g 2 )/(6g 2 ) 
_ ^--) = (1 + e6 )p(_4p 2 + g 4 - g 2 p - g 2 )/(4g 2 ) 
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+K b p(12p 2 - 3g 4 + 3q 2 p + ll<f)/(12g 2 ) 
KW+- _ W+T) = (1 + e b )p(q 2 + p- l)/(2v^V?) 

+i^(-g 2 + p- l)/(3y/2y/F) 

1{w +- + = (1 + ^{-fy^yf^) 

+K b (3p 2 + 2g 2 )/(3v / 2V / ?) 

where p = — p + g 2 ) 2 — 4g 2 . For the sake of completeness we have retained the 

contribution of the chromomagnetic interaction G b contribution, although it vanishes for 
Afc-decays. For applications of the helicity formalism to the mesonic sector one has to set 
P = in the rate formulae and replace (W?\, + v ) by W\,.,y . 



26 



APPENDIX B: FULLY INTEGRATED RATE FUNCTIONS 



In this Appendix we list our fully integrated results for those structure function 
components that determine the angular decay distribution for unpolarized A;,-decay, i.e. 
for the helicity rates f„, f f^, f f £, f J and fj L . 



1 , = (1 - - 777 - 7/f • n" ~l> ■ VVP - ><fl> ~ 7p 2 - 7 V p 2 + p 3 ) 



+16tf 6 ^— (1 - 577 - 2r/ 2 + lOp - 5rpp + p 2 ) 
9 



pV- 

(1 - 2?7 + rj 2 + p) In 



/ l-?7 + p-Vg \ 2 2 / 1 + 77 -P - VR \ 

1] 111 1 — w? — J + 77 ip ~ 1] 111 1 — — J 

1 -r] + p- VR\ 



2VP 



+-^-(p-l)ln 



( l + V-p- 
2^ 



(1-^6 



(1 + IO77 + ?7 2 — 7p — 10?7p + rfp - 7p 2 + lOrpp 2 + p 3 



^^7T^(1 — 5r7 - 2?7 2 + 10p - 5r/p + p 2 ) 



9 



+ (1 - i^)8p 2 (-2 + 3 V - rj 2 - V p) In 



1 — rj + p 



2v/p 



+K b —(-l + 2V ~ V ) In ^ — J 

+ (1 - K 6 )8r7(l - p)(l + r] - 2p + r,p + p 2 ) In 

+^(w>,(i±^) 



2^ 



7 = (-1 + 77 + 2^P - p)(l -If]- 7r] 2 + rf + 2^p- 12 V ^p - 2r, 2 ^p - Yip + 38r/p 
-7r/ 2 p + 28v/p 3 - I2r,^f - 17p 2 - 7r,p 2 + lyff + p 3 )/3 
+V(l-p) 2 ln' 77 » 



+K b 



^ (1 - 77 - 2 y/p + p) (9 - 2377 + 77 2 + r> 3 - 30^p + 2O77 ^p - 2r?Jp 



+31p + 2277p - 77/ 2 p - - 12rj^p 3 - 9p 2 - 7r,p 2 + 2jf + p 3 ) 



I677 2 



(l-p) 2 ln 



f + = (1 - K b )^^(l + IO77 + T/ 2 - 2p + IO77P + p 2 ) + K b ^f^(l + 577 + p) 
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+ (l-K b )j^-(-l + 7 1 + p)\ni '- 

+K b J 2riP ^ {l -2r] + r] 2 -p + 2rjp) In 



2 VP J 
fl-rj + p- v^R N 



3(1 -p) 



+ (1 - K b )—^—(l - p- p 2 + p* + r] + rip 2 ) In 



?Vp J 

(l +r] - p- VR 



l-p 



+K b 



32rf 



3(1 -p) 



(2 + r? - 4p + 7]p + 2p 2 ) In 



2^ 

1 + 77 - p - ^[R \ 



(67) 



f + = (1 - K b ) V VR(-3 - 3r/ + 4p - 3r 1 p - 3p 2 ) - K b ^j^(l + 5 V + p) 

+ (l-K b )p^(-3 + 4r ] -r ] 2 + 4p-4 V p-p 2 )\n(-~ V + P ~^ 



l-p 

I^ J 2VP A -l + 2r l -rj 2 + p- 2rjp) In 



2^ 
1 — 77 + p — y/R\ 



3(1 -p)" ~ ' - ' ' ' r 2 ^p 

-1 - rf + 4p - 6p 2 - 7? 2 p 2 + 4p 3 - p 4 , 



+ (l-K b )2r/(-4r/(l-p 2 



+#6 



32r/ 2 



-2 - 7] + 4p - i]p - 2p 2 ) In 



3(1 -P) 

f + = (1 - K b ) \r]VR(-3 - 3r] + 4p - 3?7p - 3p 2 



l-p 
'l + rj-p-VR s 

2 Vv 



In 



'1+77- p-Vr s 



2^ 



+2 ,(-4,(i - p 2 ) + - 1 -^ + 4 ^- 6 ^-^ 2 + V -^ 4 ) Ji+n-p-JV 



l-p 

2r/p2 -(3-4r ? + r ? 2 -4p + 4r ? p + p 2 )ln A-^ + P"^ 



2^ 



l-p 



2^ 



f + £ = ^^^(i " V ~ 2VP + P)(3 + 3/7 - 2^p - A V ^p - 2p + 3 V p - 2 ^p 3 + 3p 2 
-r/(l + 4r/ + r/ 2 - 3p - 7] 2 p + 3p 2 - Arjp 2 - p 3 ) In 



.(i-VpT 



[^(1 + Vp)(-1 + 77 + 2^p - p)(-l -77-3^ + 3VP + p + 3^P 3 ) 
+ y (P " 1)("1 - V 2 + 2p - 47/p - p 2 ) In ( (1 _^ )2 



(6* 



(69) 



(70) 



where R is related to the maximal momentum of the r in the A b rest frame and is given 
by 

R = 1 - 277 + r/ 2 - 2p - 2r/p + p 2 



As a necessary check we reproduce the total rate formula |||7]j8|,^] : 
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T = T b fa + f+ + tl + f+ + 3f+) 



r b (i-K b ) 



\/R{\ -It]- 7rf + rf -7p+ Ylr)p - 7r] 2 p - 7p 2 - 7r]p 2 + p 3 ) 



+2A7] 2 (p 2 - 1) In 



1 + 7] - p 



R 



2^ 



+ 2Ap 2 (7] 2 - l)ln 



l-r] + p-VR 
2^p 



(71) 



We mention that we have double checked all our analytic results by comparing them to a 
numerical evaluation. 
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FIGURES 




FIG. 1. Definition of the polar angles G and Qp and of the azimuthal angle <p in the decay 
A& — > X +W~{-^> in the A& rest system. We specify a z-axis which we take to be along 

px- P denotes the polarization three- vector of the Afe. 




FIG. 2. Lowest-order graph contributing to the forward scattering amplitude T(s). 
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